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Abstract. We explore some qualitative dynamics in the neighborhood 
of the 3 — dimensional two-fold symmetric singularity. We study the 
existence of an one-parameter family of regular (pseudo) periodic orbits 
of such systems near a reversible two-fold singularity. 



1. Introduction 

Theory of non-smooth dynamical systems has been developing at a very 
fast pace in recent years and it has become certainly one of the common 
frontiers between Mathematics and Physics or Engineering. Hybrid and 
switched models are being increasingly used in applications to describe a 
large variety of physical devices. Examples include mechanical systems with 
friction and backlash, electrical and electronic circuits, walking and hopping 
robots and, more recently, biological and neural systems [B-S-C] IC j IK-R-G] . 
One should observe that much work has been done in the study of the 
qualitative aspects of the phase space of discontinuous vector fields. 

We consider vector fields expressible in the form x = Z(x) where x G I 3 is 
a state vector and Z is a smooth piecewise mapping. The discontinuities are 
concentrated on a codimension-one submanifold E of M 3 . E is usually called 
the switching manifold. Orbit-solutions on E, whenever possible are defined 
according to the Filippov convention, [Fj. We point out that trajectories 
may become constrained to the switching manifold, and this behavior is 
called sliding. 

A general understanding of dynamics of generic 3D Filippov sytems was 
obstructed by the appearance of the two-fold singularity (refer to 12.21 and 
[Tip . Observe that, topologically speaking, its shape is very simple, and 
moreover it is generic in piecewise smooth systems with three or more di- 
mensions. 

In the last years we have been noticing a great interest in the study of 
the two- fold singularity (also called T-singularity) . In |J-Cj is classified the 
dynamics in a neighborhood of a T-sing ularity. In [C-B-F-Jj her occurrence 
was discussed with real models in control theory, introducing conditions for 
the study of the existence of this singularity. 
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Let us consider a discontinuous vector field Z = (X, Y) and its corre- 
sponding normalized sliding vector field Z (see Definition Q] and Remark 
[I]). Let po a point where Z presents a discontinuity. When po is the two- 
fold singularity, po is always a critical point of Z s and a fixed point of the 
return map (see definition [5]) which is a composition (p z = 7 y o 7^ of two 
involutions associated respectively to Y and X. 

We deal with perturbations Z (whose expressions will be precise in 13. ip 
of the 3D piecewise smooth vector fields Zq = (Xq,Yq). Our aim is to find 
conditions to the existence of typical closed orbits of Z (see definition [9|) . 
To do this we use some geometrical properties, like a reversibility. So, here 
we exhibited a subset of this typical singularities that is reversible. 

The paper is organized as follows: in sections [2] we deal with some pre- 
liminaries, give some definitions, and establish the notations. In section [3] 
we state the results about the reversibility. In section [3] and 0] we prove the 
existence of a family of the regular and pseudo periodic orbits, respectively. 

2. Preliminaries 

2.1. Distinguished regions of the discontinuity set. In this section 
some notations, basic definitions and elementary concepts are presented. 

Designate by X r the space of all germs of C r vector fields on R 3 at 
endowed with the C r -topology with r > 1 and large enough for our purposes. 
Call Q r the space of all germs of vector fields Z in (R 3 ,0) such that 

7ln s _ / X(q) t for h(q) > 0, 
W ~ \ Y(q), for %)<0, 

where h is (a germ of) a smooth function h : (]R 3 ,0) — > (R, 0) having 6 R 
as its regular value. Let £ = h~ 1 (0). We assume that G £. 

The above vector field is denoted by Z = (X, Y). We endow £l r = X r x X r 
with the product topology. 

To define the orbit solutions of Z on the switching surface £ we follow 
a pragmatic approach. In a well characterized open set O of S (described 
below) the solution of Z through a point p G O obeys the Filippov rules (see 
[F]) and on £ — O we accept it to be multivalued. As we are dealing with 
systems derived from ordinary differential equations the non-uniqueness of 
solutions is allowed. We just must take into account all the leaves of the 
foliation in R 3 generated by the orbits of Z (and also the orbits of X and 
Y) passing through or exiting from or converging to a point p G £. 

For each X G X r we define the smooth function Xh : R 3 — > R given by 
Xh = X.\7h where . is the canonical scalar product in R 3 . 
In what follows we use the Filippov convention (refer to [F]). We first 
distinguish the following regions on £ : 

► Sewing Region: SwR = {p G £; (Xh)(p)(Yh)(p) > 0}. When 
convenient we denote SwR + = {p G £; (Xh)(p) > 0, (Yh)(p) > 0} 
and SwR~ = {p G £; {Xh)(p) < 0, (Yh)(p) < 0}. In general a point 
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in the phase space which moves on an orbit of Z and reaches a point 
in SwR, crosses S. 

► Escaping Region: EscR = {p £ E;(Xh)(p) > 0,(Yh)(p) < 0}. 
In this case any orbit which meets EscR remains tangent to £ for 
negative times. 

► Sliding Region: SIR = {p £ E; (Xh)(p) < 0,(Yh)(p) > 0}. On 
SIR the flow slides on S; the flow follows a well defined vector field 
Z s called the sliding vector field (see Definition [T]) . 

Generically, the set O = SIRU EscRU SwR is open and dense in S. Observe 
that for any p £ 0, we have X(p) ^ and Y(p) ^ 0. 

Definition 1. Let Z = (X, Y) £ fl r . The sliding vector field Z s associated 
to Z is a linear convex combination of X and Y tangent to S ; that is, 

Remark 1. Observe that EscR for Z represent SIR for —Z. We can define 
the escaping vector field on £ by —{—Z) s . This vector field is called the 
sliding vector field independently of whether it is defined in the sliding or 
escaping region. For p G SIR U EscR the local orbit of p is ruled by this 
vector field. Therefore all this orbit is contained in SIR U EscR, and the 
future orbit of the sliding vector field coincides with the future orbit of the 
normalized sliding vector field Z s = (Y.VhX — X.VhY). 

Observe that Z s and Z s are orbitally equivalent on SIR (resp. on EscR). 
Moreover Z s can be C r -extended beyond the boundary of SIR, EscR. For 
technical reasons we consider the future orbit of Z through a point p £ 
SIR U EscR given by the orbit of Z s . 

Notation: In all what follows we consider the map h : (x, y, z) i— > z. So, 
the expression of Z s is: 

Z s = {X l Y* - Y l X z , X 2 Y 3 - Y 2 X 3 ), 
where X = (X\X 2 ,X 3 ) and Y = (Y\Y 2 ,Y 3 ). 

Definition 2. We say that is a two-fold singularity of Z = (X, Y) G £T r 
if Xh(0) = Yh(0) = and X 2 h(0) ^ 0,Y 2 h(0) ^ 0, where X 2 h(0) = 
X(Xh)(0). 

Definition 3. If p £ SIR U EscR and X(p),Y(p) are linearly dependent 
then p is a critical point of Z s . In this case p is called a pseudo equilibrium 
ofZ. 

The curves of tangential singularities or the S— singularity of X in S are 
given by Sx = {p G S; Xh{p) = 0}. 

We denote by Q, F C £l r the set of non-smooth vector fields such that the 
origin is a two- fold singularity. First of all, observe that: 
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1- The trajectories of both X and Y through have a quadratic contact 
with £ (at 0). 

2- Generically, Sx and Sy are transverse at 0. This case was first 
studied in [Tlj . 

3- If is a two-fold singularity of Z then Z s can be C r -extended to a 
full neighborhood of G S, Z s (0) = and tp z (0) = 0. 

2.2. Two- fold singularity. The following construction is given in |T2j . 
Let Z = (X, Y) G O r such that is a fold point for both X and Y. Applying 
the Implicit Function Theorem, for each p G (£, 0) there exists a unique t(p) 
such that the orbit-solution 1 1— > <px(t,p) of X through p meets S at a point 

p = 4>x{t{p),p)- 

We then define the smooth mapping 7x : (M 2 , 0) — > (M 2 , 0) by 7x(p) = V- 
This map is a C T '-diffeomorphism and satisfies: 7^ = Id. Analogously, we 
define the smooth map associated to Y: 7y : (]R 2 , 0) — > (M 2 , 0) which satisfies 
7y = Id. We define now the first return map associated to Z = (X,Y): 

Definition 4. The first return map (fz '■ (£, 0) — > (S, 0) is defined by the 
composition 

ip z = 7 y O 7 X 

H^e denote by Lz(., ■) the linear part of ' <pz- 

(pz is a C r -diffeomorphism preserving area: det(Dipz)(0, 0) = 1. So the 
eigenvalues of Dip(0, 0) are /3 and /3" 1 with: 

(a) saddle type: /3 G R, /3 / 0; 

(6) elliptical type: j3 = e ie with g]0, 7r[. In this case, Lz(-,-) is a 
rotation. 

2.3. Partition of fi^. Consider the subset of 0^: 

Regular two-fold: Let Qq {5) be the set of all Z = (X, Y) G n F 
such that the contact between Sx and Sy at is transverse, the 
eigenvectors of DZ s (0) are transverse to Sx and Sy at and is a 
hyperbolic critical point for Z . 

Remark 2. From [S-T] we derive that £Iq(6) is a codimension zero subma- 
nifold of ft. 

In £Iq(5) we distinguish the following subsets, corresponding to three 
cases: 

Elliptic case: Sift (6.1) = {Z G Qfi(8); X 2 h(0) < and Y 2 h(0) > 0}. 
We have two invisible tangencies (invisible two-fold or a T- singularity). See 
Figure [TJ 

Parabolic case: n£(<5.2) = G 0^(<5); A 2 /i(0) > 0,y 2 /i(0) > or 
X 2 h(0) < 0,Y 2 h(0) < 0} (visible fold- invisible fold); 

Hyperbolic case: tl$(6.S) = {Z G Q% (5); X 2 h(0) > 0,Y 2 h(0) < 0} 
(visible two- fold). 
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s x S x 
Elliptic case Hyperbolic case 



Figure 1. Two-fold singularity. 

Recently, many works tried to understand the dynamics of the T-singularity, 
|J-Q lJl [C-B-F- J] . See |T2j for further references and related topics. 

Definition 5. We say that a C°° -diffeomorphism £ : HL 3 — > M 3 , is an invo- 
lution if £ o £ = Id. We put Fix(£) = {p 6 K 4 | £(p) = p}. 

Definition 6. We say that a non-smooth vector field Z G Q r with discon- 
tinuity manifold X is ^-reversible if there exists an involution £ : M 3 — > R 3 
such that: 

1- Fix(Q C S 

2- VpG'S, toZ{p) = -Z(i{p)). 

Definition 7. Let ^ e ^ e se ^ °f dements Z G Qq(o~) such that: 

1- Z is ^-reversible; 

2- Z is a generic invisible two-fold at 0; 

Definition 8. We say that Z G O r is simple at if X 2 h(0) / ±XYh(0). 

We say that X G 3L r is semi-linear if X l (x, y, z), X 2 (x, y, z) possess only 
the 0-jet (degree zero) and X 3 (x,y, z) possess only the 1-jet (degree one). 
Put: 

SL(n% (5)) = {Z = (X, Y) G (5);X and Y are semi-linear}. 
Consider the set of non smooth dynamics vector fields: 

QP(8) = {Z G SL(fi£ (5)); X 2 h(0) < 0, (X(Yh))(0) < 0, 

(X(Yh))(0)(Y(Xh))(0) = {X 2 h)(0)(Y 2 h)(0),Y 2 h(0) > 0}. 

2.4. Orbits. By convention, if p G SIR the future orbit of Z through p is 
given by the trajectory of the sliding vector field Z s through p. 

Definition 9. Let Z = (X,Y) G O r and x G (M 3 ,0). 
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1- The orbit t i— > <j)z(t,xo) is a regular periodic orbit if it is closed 
and composed by segments of orbits of X, Y and Z s keeping the 
orientation. See Figure® 

2- We call pseudo periodic regular orbit a trajectory of Z which is 
closed, composed by segments of the orbits of X and Y, with non 
preserved orientation (see Figure^). 




Regular Pseudo 
Figure 2. Periodic orbits. 

2.5. Main results. Consider the set of straights line passing for the origin: 

Ca = {(x,ax,0) £ S; with a £ R}. 
We obtain the result about the existence of pseudo periodic orbits and 
periodic orbits for Z £ ^o(^) : 

Theorem 10. Let Z £ fig (5) and p £ C a . Then: 

(a) If Z is simple at then Z does not have k-periodic orbits for all 
k £ N. 

(b) If Z is non simple at and a = 1 then Z has one-periodic orbit 
passing through p. In addition, for a / 1 does not exist k-periodic 
orbits for all k £ N. 

(c) If Z is non simple at and a = — 1 then Z has one-pseudo periodic 
orbit passing through p. In addition, for a ^ — 1 does not exist 
k-pseudo periodic orbits for all k £ N. 

(d) is a hyperbolic equilibrium saddle (elliptical) point of ipz provided 
that \XYh(0)\ > \X 2 h(0)\, (\XYh(0)\ < \X 2 h(0)\), see Figured 

In Theorem 1101 we exhibit conditions on Z £ Og(<5) for the existence of 
one-parameter families of periodic and pseudo periodic orbits for Z. 

We also study the existence of pseudo periodic orbits and regular periodic 
orbits for Z £ SL(Qq (5)). Observe that hypothesis in Theorem [y] (a) 
complete the result obtained in Theorem I1UI (c) for Z £ SL(Qq (5)). In the 
item (b) we get a more general family of regular periodic orbits, generalizing 
the family obtained in Theorem 1101 (b), for Z £ SL(Qq (5)). There results 
are summarized in: 
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Theorem 11. Consider Z = (X,Y) G SL{n£(5)),p G C a and X 2 h(0) < 
0,y 2 /i(0) > 0. Then 

(a) IfZ is simple at the origin, XYh(0) = -YXh(0), X 2 h(0) = -Y 2 h(0) 
and a = — 1 then Z has a 1-pseudo periodic orbit passing through p. 
In addition, for a/-l there does not exist k-pseudo periodic orbits 
for all fceN. 

(b) If Z G Q,q(5), a = X* and p G SwR then p is either a fixed point for 
(pi or ip 2 , where A* = xYh(o) ■ 

3. Reversibility 

In this section we will explore some symmetry properties of Z G 
Our objective is to study the existence of periodic orbits of arbitrary period 
k. 

3.1. Normal forms for the reversible two- fold singularity. Consider 
Zq G Qq(5). Throughout this section, we fix a coordinates system (x,y,z) 
such that 

€(x,y,z) = (y,x, -z). 

Observe that Fix(£) = {(x, y, z); x = y, z = 0}. 

Let / : M 3 — > K be a non-zero polynomial in (x,y,z). We call dr(f) the 
total degree of /, that is the maximum sum of exponents of the monomials 
of/. 

Consider the polynomials: f x a : M 3 — > K where a G {x, y} and i = 1,2, 3, 
with the following specifications: f a ' have null linear part and dx(f^) > 1- 
Let 

F a (x,y,z) = (fl(x,y,z),fl(x,y,z),fl(x,y,z)). 

In what follows we exhibit the topological normal form for the two-fold 
singularity: 

Proposition 12. (Normal Forms) If Z G then Z is C Q — equivalent 

to Z = (X, Y) where 

X(x,y,z) = (C x ,Cxy,x) + F x (x,y,z), 
Y{x,y,z) = {C Y x,C Y ,y) + F y (x,y,z), 

with C x = X 2 h(0),C Y = Y 2 h(0),C X Y = X(Yh)(0), C Y x = Y(Xh)(0), 
X 2 h(0) / and Y 2 h(0) / 0. 

Proof. The origin belongs to the frontier of the sliding region. Consider 
a local coordinates system such that Sx = {(x, y, z); x = z = 0} and S Y = 
{(x,y,z);y = z = 0}. With these settings the local normal forms of X and 
Y are: 
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X(x,y,z) = (C x ,C XY ,x) + F x {x,y,z) 

Y{x,y,z) = (C YX ,C Y ,y) + F y (x,y,z). 

Observe that X 2 h(0) = C x ,Y 2 h(0) = C y ,C X y = X(Yh)(0) and C YX = 
Y(Xh)(0). 

m 

Consider the regions in E: 

SIR = {(x,y,0);x <0,y>0}, 
(1) SwR ={(x,y,0);x.y > 0}, 

EscR = {(x,y,0);x > 0,y < 0}. 

3.2. Dynamic of (pz when Z is £— reversible. The next Lemma exhibit 
the subset of that is ^-reversible. 

Lemma 13. Let Z G Sl$(5). If XYh(0) = -YXh(0), X 2 h(0) = -Y 2 h(0) 
then Z is ^-reversible. 

Proof. Straight forward computations. 

■ 

The proof of Theorem [10] item (d) follows by: 

Lemma 14. Let Z = (X,Y) G fi{j(<y). Then: 

(i) is a hyperbolic equilibrium saddle point oftpz provided that \XYh(0)\ 

> \X 2 h(0)\, see Figure^ 
(ii) is a hyperbolic equilibrium elliptical point of (pz provided that 
\XYh(0)\ < \X 2 h(0)\. 

Proof. Consider the local normal form of Z = (X, Y) given in the 
Proposition Q2J By the Lemma US we have: XYh(0) = C XY ,YXh(0) = 
C YX ,X 2 h{0) = C x and Y 2 h(0) = C Y . 

The regions in E are given in ([TJ. Let po = (xo,yo, zq). The flows of the 
vector fields X and Y are: 

X (Po) = {x + C x t,y + C XY t,z + x t + ±C x t 2 ) +0(t 2 ,t 2 ,t 3 ) 

(2) 

<j) Y (p ) = (s + C YX t,y + C Y t, z + y t + \C Y t 2 ) + 0(t 2 ,t 2 ,t 3 ). 
If po = (so, 2/0,0) G E, taking 

, , x 2x !»/ \ -4C X yx + 2C x y 

*i(Po) = -7^- and t 2 {po) = 

<^x C x 

we obtain <p x (po) = (xi,yi,0) G E and 4> Y {%\ , 2/1 , 0) G S. We define the 
return region (the shaded region in the Figure [3]) by: 

R = [{(x,y,0) G E; x > 0, y < 2C , X yC^ 1 x} n Su;.R]. 
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The return region Rq represent the points p = (x,y,0) S X such that 
U(p),i = 1,2 are positive. The first return map is expressed by: 




Elliptical type Saddle type 

Figure 3. Dynamics of ipz- 



(3) Mp) = # o #(p) = (( 4Ci ^ Ci ) * - 2 -^y, ^ - v) ■ 

The eigenvalues of ipz are A± = — 1 H — — 1 — — — ^ . 

Besides, the origin is a critical point of the saddle (respectively elliptical) 
type of if Z if \C X y\ > I C x | (respectively \C X y\ <\C x \). 



3.3. Existence of k— periodic orbits for Z. We investigate the existence 
of periodic orbits of Z £ 0q(<5) of period k arbitrary. We prove the items 
(a) , (6) and (c) of Theorem [TO) 

Proof. By ©, the equation ipz(x,y,0) = (x,y,0) is expressed by: 

~~ C x + ^C^y 2Cxy 

c\ X ~~c x ~ y =x 

2C X y 

— x-y =y. 

Gjsf 

The solutions are C X y = ±C X . So, the periodic orbits pass through the 
pair of straight lines {(x,y);x = =ty}. Observe that the straight line 

C_i = {(x,y,0);x = -y} 

is contained in SIR U EscR. So we get a one-parameter family of pseudo 
periodic orbits. 

On the straight line 

Ci = {(x, y, 0); x = y} = Fix (£) 

we obtain a one-parameter family of periodic orbits. The expression of ip^ 
is given by: 
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ip n (x, ax) = <Pz(x } ax) =( 



-c x +ac 2 xy 2aCxY 



ci 



C x 

n-l 



X, 



-C x +iC XY _ 2aC XY 

Solving ip n ( ) = (x, ax), we obtain the system 
-C X + 4C 2 XY 2aC XY " 



2Cxy 
Cx 



a 



x). 



r 2 



C x 



-C\ + AC\ Y _ _ 2aC X Y 
C x 



-, n-l 



r 2 



2Cxy 
Cx 



a 



a. 



We obtain directly that Cx = 2 °+q/ . Replacing this value in the first 
equation, we obtain a~ 2n = 1, ie, a = ±1. So the conclusion of the Theorem 
[10] is straight forward. 



4. PSEUDO PERIODIC ORBITS 

4.1. Existence of A;-pseudo periodic orbits for Z. A necessary condi- 
tion for the existence of pseudo periodic orbits is 

(fi(x,y) = <p 2 (x,y) 

where <pi(x,y) = (<$ o cf> x ) and (p 2 (x,y) = (4> x $y )( x ' V-> °)> witn ^ °> 
for i = 1, . . . , 4. See Figure HI 

The regions on £ are given in (pQ). We define the subset of EscR C S: 

i?5 - I fx „ 0) G EscR- 2Y (Xh)(0) X 2 h(0) 
RS - | (x, y, 0) G EscR, y2/i(Q) 2/ < * < ^YT^of 

Lemma 15. lei Z G SL(fl$(S)). If X 2 h(0) < 0,Y 2 h(0) > 0,XYh(0) = 
-YXh(0),X 2 h(0) = -Y 2 h(0) and XYh(0)YXh(0) = X 2 h(0)Y 2 h(0) then 
RS = H>. 

Proof. In fact, 

2YXh(0) X 2 h(0) 



-y < 



y ^ 4X 2 h(0)Y 2 h(0) > X 2 h(0)Y 2 h(0). 



Y 2 h(0) u 2XYh(oy 

But this is a contradiction since X 2 h(0) < and Y 2 h(0) > 0, by hypothesis. 

■ 

We recall the notations: 



C a = {(x,ax,0) G £; with a G 
Now we conclude the proof of Theorem [IT] 
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Proof of Theorem [TT] Consider the regions in E given by ([I]). The flow 
given in (|2j). So, the expressions of (pi(x,y) and ip2(x,y) are: 



where h(x,y) = -^,t 2 (x,y) = -^(-^i + y),t 3 (x,y) = --^y and 
t A (x,y) = - jjfix - ^y). As X 2 h(0) < and Y 2 h(0) > we get h > 
and *3 > 0, respectively. 

To prove item (a) we need to solve <pi(p) = <P2(p) with p = (x, y, 0) G E. 
Initially we suppose that p G EscR. From @ and solving (pi(x, y) = f2(x, y) 
we obtain: 

(i) (x(Yh))(o)(Y(xh))(o) =jxy(o)(Y 2 h)(oy, 

{%%) (x,y) G C\*, where A* = xYh(o) ■ 

As C\* C EscR we obtain that XYh(0) > 0. Replacing this condition in 
(i) we get YXh(0) < 0. 

In this way, we consider the diffeomorphisms (pi and q>2 restricted to 
RS C E. By Lemma 13 i?5 = 0. 

Therefore, for Z G SL(Uq(6)) there does not exist pseudo periodic passing 
by p G EscR. 

If p G Swi? follows by item (c) of Theorem 1101 the result. Observe that in 
this case, some values of the time ti, i = 1, . . . , 4 are equal to 0. 



Figure 4. In (a) is represented a pseudo periodic orbit. In 
(6) is represented one of the periodic orbits of Z G SL(Q). 

Let us prove item (6). For Z G Qq(5), the conditions (i), (ii) and = ip 2 
are satisfied since C\* C SwR. The dynamics of Z G is illustrated in 

Figure 0} 

Claim: "If p G Cj. C Su/-R + and Z G SIq{5) (respectively p G C A * C 
SwR~ and Z G then = p (respectively </?2(p) = p)- In other 

words, on C\* there is a family of periodic orbits" . 

In fact, observe that if Z G Qq(5),p G C\* C 5wi? then Z satisfies 
the conditions (i), (ii) and t{ > 0, for i = 1 ... ,4. That is, (fi(p) = p (p G 
C\*nSwR + ) or (/?2(p) =p(p£ C\*nSwR~). Besides, for Z G ^((J) we have 




(4) 




A" 



X 
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C x * C (SwR+USwR-). If p G [C x * nSwR-} then ^(p) = $?o0%(p) = -p 
and y) = 0. 

■ 
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